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We investigate the non-Gaussianities of inflation driven by a single scalar field coupling non- 
minimally to the Einstein Gravity. We assume that the form of the scalar field is very general with 
an arbitrary sound speed. For convenience to study, we take the subclass that the non-minimal 
coupling term is linear to the Ricci scalar R. We define a parameter /i = eh/eg where and 
eg are two kinds of slow-roll parameters, and obtain the dependence of the shape of the 3-point 
correlation function on /i. We also show the estimator Fnl in the equilateral limit. Finally, based 
on numerical calculations, we present the non-Gaussianities of non- minimal coupling chaotic inflation 
as an explicit example. 
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I \ 

! I. INTRODUCTION 
(N . 

The inflation Theory is one of the most successful theories of modern cosmology. Having a period of very rapidly 

1 I ' accelerating expansion, it can not only solve many theoretical problems in cosmology, such as flatness, horizon, 

• monopole and so on, but also gives the right amount of primordial fluctuations with nearly scale-invariant power 
' spectrum, which fits the data very well in structure formation [iHl]. 

. There are many ways to construct inflation models, one of which is to introduce a scalar field called "inflaton" (f) (see 
Moreover, one may expect that inflaton could have non- minimal coupling to Ricci scalar R. The most usual 
I""!. coupling form is R(j>^ which was initially studied for new inflation scenario [3] and chaotic inflation scenario {sj. Later 
'"pi ■ on, various models have been taken on with deeply and wildly investigations. With a non-minimal coupling term, 
' ' inflation can be easily obtained and an attractor solution is also available ^]. Perturbations based on non-minimal 
coupling inflation are discussed in Q , where the coupling term may give rise to corrections on power spectrum which 
(— I ] can be used to fit the data or constrain the parameters. Non-minimal couplings can be extended to multifields, see 
'— Q, or kinetic term coupling [9|. The constraints from observational data were also performed, e.g. in [T^], where the 
■ authors claimed that for non-minimal coupling chaotic inflation models, a tiny tensor to scalar ratio will be obtained. 
^ , Other applications of non-minimal coupling inflation include the realization of warm inflation |1I|] and the avoidance 



J — ■ of the so called "77" problem [12| in the framework of string theory [13|. One can also see [l4| and also [15| for 
^\ ' comprehensive reviews of non-minimal coupling theories. 

\^ The non-Gaussianity of the primordial perturbation has been widely acknowledged to be an important probe in 

■ the early universe |16l - l22| . Experimentally, more and more accurate data allow us to study the non-linear properties 

■ of the fluctuation in Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) [23l - l25l |: Theoretically, 
T-H the redundance of inflation models requires more information than those of linear perturbations only to have them 

distinguished. The non-Gaussianity of the fluctuations was first considered in [26|, and it was further shown in [l^] 
^""^ ' that the canonical single field slow roll inflation can only give rise to negligible amount of non-Gaussianity. To get large 
^ ' non-Gaussianity peogle need to find new inflation models, an INCOMPLETE list and references of which_^include: 
' multi-field models k-inflation [2^, DBI-type inflation [sTl, curvaton scenario [s^, ghost inflation [32], warm 
rN \ inflation j33|], non-Bunch-Davies vacuum scenario l34l . bounce scenario [IH], island cosmology [36j, loop correction 
. [131, non-commutativity Issl . string gas scenar io |40|| . cosmic string [s^], "end-in- inflation" scenario jij, Ekpyrotic 
■ ■ ■ ' scenario [421], vector fleld ^3"], Horava theories j44l] and so on and so forth [Glj. 

In this note, we investigate the non-Gaussianity of inflation driven by a general single field P{X, (p) coupling non- 
minimally to the Einstein Gravity. Some specific examples of non-Gaussianities of non-minimal coupled field has 
been studied in, e.g., [i^ and non-Gaussianity generated by modified gravity is expected to have effects that can be 
tested by CMB anisotropies [13]. By taking a subclass of linear coupling, we calculated various shapes depending 
on the ratio between two slow roll parameters et and eg, which describe the evolution of cosmic expansion and the 
non-minimal correction, respectively. The power spectrum will deviate from scale- invariance due to the existence of 
non-minimal coupling 15], and the shape of the 3-point correlation function are correspondingly affected. In this 
paper we find that for different (red or blue) tilt of the power spectrum, the shape will include different parts which 
will obtain different amplitude of non-Gaussianities. However, since we have only calculated up to leading order in 
the slow-roll parameter, this conclusion has not been so unambiguous yet. Nevertheless, If it can be verified after a 
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complete consideration to all the orders, one can find the relations between 2- and 3-point correlation functions which 
can be used to constrain non-minimal coupling models. This will be one of our future works. 

This paper is organized as follows: Sec. II briefly reviews the preliminaries and basic equations of the general non- 
minimal coupling single field inflation. We study the non-Gaussianities of the general non-scalar field with linear 
coupling in Sec. Ill, which is the main part of the paper. We first study the perturbed action of the system up 
to 3rd order, and obtained the mode solution at the quadratic level. After that, we calculate various shapes of the 
3-point correlation functions using the mode solution. We also study their equilateral limit and the relation with 
slow-roll parameters at their leading order. In the last part of this section, we present the non-Gaussianities of non- 
minimal coupling chaotic inflation as an explicit example using numerical calculations. Sec. IV is the conclusion and 
discussions. 



II. PRELIMINARY 



To begin with, let's consider the most general action of a single scalar field with non-minimal coupling: 

S^IJ dtd^xV^ifiR, 0) + 2P{X, </>)] , (1) 

where X = ~\g^'^ d^4>d,j(t) is the kinetic term and the metric g^y — diag[—l,a^{t),a?{t)^a?{t)\ with a{t) the scale 
factor of the universe. For the background evolution, one can vary the action ([1]) with respect to the field and the 
metric to get the equation of motion for </>: 



and the Einstein Equations: 



where 



2P^ + 2{PxxV''X + Px0V^0)V^0 + 2Pxn0 = , (2) 
= T^^^ , (3) 



Spi. = Ofngf,^ - \7^\/^fji + fjiRf^t, - ]^fg^u , (4) 

T^^) = PxV^0V,0 + Pg,,, , (5) 

with being the covariant derivative with respect to the metric (7^1, and □ = V^V. 
The evolution of the universe can be described by the slow roll parameter: 

where H = d/a is the Hubble parameter and dot means the time derivative. In the inflation case, we require that e/j 
be small. Moreover, we can define two more parameters following p7| : 

E = XPx + 2X^Pxx , (7) 
A = X^Pxx + ^X^'Pxxx , (8) 

which will be used in later parts of the paper. 

To study non-Gaussianities, we adopt the usual convention of using the Arnowitt-Deser-Misner (ADM) metric [i^ 
as follows: 

ds^ = -N^dt^ + h,,,j(dx' + N'dt){dx^ + N^dt) , (9) 

where N{x) and Ni{x) are the lapse function and the shift vector, respectively. It is useful to decompose the action 
^ into the 3+1 form, saying: 

S ^ \ I dtd^xy/hN{f[R^^'> + KijK'^ - + -^dt{VhK) '^-=d,{VhKN' + Vhh'^ dj N) , cji] 

2 J Ny/h N\/h 

+2P(A,0)}. (10) 
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The extrinsic tensor Kij is defined as: 

K^J = ^{k, - V.iV, - V.iVO (11) 

where is the covariant derivative with respect to the metric hij and its indices can be raised and lowered by h^j. 
The contraction K = K]. The three-dimensional Ricci scalar R^'^^ is computed from the metric hij. From this action, 
we are able to obtain the equations of motion for N and Ni (constraint equations) as: 

-^djiVhh'^djR) +2P- 2Pxv^N-^ = , (12) 

VjifRK"^) - K'^VjUrK) - h'^VjiN-^dtfB) + h''Vj{N-^N'difR) + h'^d.fRK - a'^ N'' Pxvd'cf> = , (13) 

respectively, where v = N^di<f). Moreover, a gauge choice is needed to eliminate the redundant degrees of freedom. 
Here we choose the uniform density (comoving) gauge, where the perturbations of the scalar field and metric take the 
following form: 

50 = , hij = a^e^^S.j . (14) 

We also need to expand the constraining variables N and Ni. It is only needed to expand them up to the {n — 2)-th 
order when we calculate n-th order perturbation 27]. So here we expand them to 1st order as follows: 

N ^l + a , N,^N,+d,^p, d'N, = , (15) 

where a, Ni and "tp are in the first order of e. 



III. THE NON-GAUSSIANITY CALCULATION 



Since it is very complicated to consider the full types of non-minimal coupling theory, in this paper we will only 
pick up a subclass where the Ricci scalar is coupled linearly to the scalar field, saying f(R, cj)) — RfR{(j)). This type 
of non-minimal coupling is often referred to as Scalar- Tensor Theory or linear coupling [H^l- For the non-linear 
coupling case, the constraint equations are listed in Appendix A and the solutions will be postponed to the future 
studies. 



A. Up to 3rd Order Action 

Considering all the equations from to (|15l) in the linear coupling case, one can get: 

-la-^fRod^N' + a-H2fR,oH + fR„)d'a - 2a"Vi?o9X = , (16) 
4al] + l2HfRoC - UH^fRoa - Aa-^fRoHd^^j - Aa-^fRod^ - {l2aH - 6C + 2a-^d^^P)fR.o = , (17) 

where /rq is the background value of /r and for the current case /rq — Jr. The equations above give the specific 
solution at first order in C,: 

TV' = , " - ^^^'^ 



2fRoH + fRQ 
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^ -2fR0 +^fR0^ g-2^ (^g) 

2fRoH + 2fRo {2fRaH + fR^f 

We define = ^ \n{fRoa^), so it will be followed that = J^qU and 6 ^ H + ^f^- Note that when Jrq 1, the 

system will return to the minimal coupling case and 9 coincides with the Hubble parameter. It is also convenient to 
rewrite a and "0 in terms oi 6 as: 

a^i, V = -7+X, X^[5a^a-ff + a'0-'e-'<^J:]d-H- (19) 
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Hereafter, we will use this parameter for convenience throughout the paper. 

Expanding the action (flUl) to 3rd order of C and substituting Eq. into the action, one can get the expanded 
action in each order: 

S'o = iy"dic;W(/o + 2Po) , (20) 
5i = , (21) 

S,^J dtd'xae''{3{^ -2^ + 1)C^ + a^e^^^E^^ + a''{^ 1 + l^m)'} , (22) 



612 

+a-4(9V)'} + 



/■2 /-a 
did3a;a3{6E4-C - 2(S + 2A)4-} 



612 



6»3 



(23) 



This process is very straightforward, but rather tedious. The physical meaning of each equation is easy understand- 
ing: The 0th order expansion ([20]) is just the background part of the original action ([T]) where subscript '0' denotes 
the background value; the 1st order expansion (1211) is just the background equation of motion. The 2nd and 3rd order 
expansions are only deviate from the GR case due to the difference of the parameter from , which will coincide 
when /flo 1- However, as will be seen below, this deviation causes very different results of non-Gaussianity of our 
case from that of GR. 



B. Quadratic Part: Mode Solution 



First of all, let's consider the solution of the 2nd order action (|22|) . This is the most important and kernel step in 
the calculation of the bispectrum which will be performed later. The 2nd order action can be written as: 



dtd^xae^^iii 



i/2 



+ 1)^2 + a2e-2«s£ + a- 



1 



= / drd 



H 



(24) 



where in the second step we have used conformal time t = J and transformed the variable Ci^j into its Fourier 
form, namely: 



It is convenient to define another variable Vk as Vk = zu-^{t)^ where 



3e2«(^ - 1)2 + ^ 



(25) 



(26) 



to let the equation of motion for Vk be in a canonical form, which is: 



Here the effective sound speed squared is defined as 



H 



HO -9)/ 



3{H- 



(27) 

and thus one can have 
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In the inflation period, the Hubble parameter changes slowly so as to have enough fast expansion. Here in order to 
solve the equation above, we need to introduce another parameter: 

ee^-l, (28) 

which describes the variation of Jro with respect to time. In case that Ch (see Eq. ([6])) and ee are both small and 
^ ei(i = h, 9), we can take the leading order of €h and eg so that Eq. ([27]) becomes: 

4 + (^'^' - - , (29) 

where n = e^/eg is the ratio of the two slow roll parameters. The solution of Eq. (I^Hl) can be presented in form of 
Hankal function: 



v^{t) = C7y^7?±(^,+ i)(c.fc|T|) , (30) 

where C is an undetermined constant denoting the amplitude of the solution. In deriving this, the approximation of 
slow- varying sound speed Cs << 1 is also taken. 

The solution (pO| can be splitted into two limits corresponding to the subhubble and superhubble regions respec- 
tively. In the superhubble region, we can take the limit to be: 



r(/i + f) r(-/i+2) 



and from the relation Vk = zu-^{t) one have: 



^ACsk\T\r+' + {c,k\T\)-^ , (31) 



.^lee + i-f I r(M+2) n-^^+2) 



where /fl,o(— 1) denotes the value of fm at r = 1 and //jo can be parameterized as fm ~ /i?o(— 1)|tP'-^^^''. In deriving 
these equations we also used the approximation t = —\/{aH) -\- O(e^). The solution contains a constant mode and 
a decaying mode, the latter of which is irrelevant and should be discarded. In the subhubble region, we can take the 
limit as: 

^ C^le-'^l^le'^ , (33) 

where we also discarded the +{^+^) branch. On the other hand, one can use WKB method to calculate the subhubble 
solution of Eq. ([2^ . which is: 

Comparing Eqs. (p3|) and ([34|) at r — oo, one can determine the coefficient C as: 



With this in hand, we can have the exact solution of u-^ in both superhubble and subhubble limits. Substituting this 
back to ([5^ . we get the final form of the superhubble solution: 

- ^ • (36) 

2/|o(-l)cr'fc^|e. + M - l|^r(-^ + i) V ^sk 
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It is a time-independent mode and thus can be applied to far future where r = 0. From this we can also obtain the 
power spectrum of C, which is: 



(37) 



and the spectrum index is = '^dlnk + 1 ~ 2(1 — /i) + 1. One can see from this that the power gets a red spectrum 
{uq < 1) when fi > 1 while a blue spectrum (n^ > 1) will be obtained at /i < 1. Moreover, the constraints that the 
primordial spectrum must be nearly scale invariant requires that \fj, — 1| ^ 
Furthermore, using Eq. ([33)) . the subhubble solution can be solved as: 



2v^c||A:|i|e^ + ^-l|2 



2^/2c|k2fl^\eh+^l-lY■ 



'^^^ + M)|Tr-i - 2»c,A:|rr] , (38) 



2V2c|fci/|o(-l)|e;, + M-l|^ 



and 



* v^/|„(-l)|« + „-l|' 

where we keep only the leading order terms in terms of e. The results above will be useful for our analysis of 
non-Gaussianity in the next paragraph. 

C. Cubic Part: Non-Gaussianities 

According to "in-in" formalism [5ll |. the 3-point correlation function is characterized in the interaction picture as: 

< ic(r,fct)c(r,fc^)c(T,fc;^)i >=-zr / dt' <\[atX)c{tX2)c{tX3),Hf„t{t')]\> , (40) 



where Hf^^ is the 3rd order interaction Hamiltonian and T is the time-ordering operator. From the 3rd order action 
(|23p. we can write down the 3rd order Hamiltonian as: 

= / dtd''xae^'{[a^e-^\j: + 2A) + 3{e - Hf]t - 3[a^e-^'j: + 3(9 - H)^]^ - a-\^ - 1 + -^)C(5C)' 

[3^2(1 _ + a^e-^e-^'nCidCf + 2a-^[3a'{l - + a^Q-^e-^'ntdCdx} , (41) 
6 9 

or, if changed to momentum space, 

d^pid^P2d^P3 . . -^,2en„2„-2e/^ , o^^ , ita rr^2l C(^, Pi)C(i, P2)C(i, pl) 



Hfnt = J ^ i^-mpi + pl + p|)ae^n[a^e-^^(S] + 2A) + 3(^ - 



(2^)9 ^ ' ■ -o. V ... V .J ^3 

or 2 -2ev. , Q//i rrN2i C(i,Pl)C(i,P2)C(*jP3) -2/^ i , ^ x/^ -^\^U -^\/■|-^ 



^a-^9-'[3a\l - + a4^-2e-2eq(^ . ^)C(t,^Jt)C(t,^)C(t,p|) 

+2a-4[3a2(l-^)2 + a40-2g-2eq(^.^)^(^^^)^(^^^)^(^^^)} ^ (42) 
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There are five terms of 3rd order, each containing a long prefactor. Using Eq. (|40p . we can calculate their contributions 
to non-Gaussianities. Neglecting the detailed calculating process, we only give the final results of the contributions 
to non-Gaussianities from each term as follows. 
The contribution from i^'^: 

+C.C., (43) 
The contribution from CC^: 



+C.C., (44) 
The contribution from C(9C)^: 

+2perms.] + C.C., (45) 
The contribution from ({dC)"^: 



l)Vi^o( l) ^_,(o)^_,(o)^_,(o) r |r|i-^^dr(27r)3<5=^(^fct)[(fc^.fc^)-^<^-.(r)u:^-.(r)^:g(r) + 2perms] 



dr -fei 

The contribution from (^dCdx- 

-'*^I^f^£fT)"^(«K-(oK^(o)/_J-I^ 



-2i^u^(0)wg(0)u,^(0) £ ^dT(27r)3j3(^fct)[(A;^-:fc^):Jii:^(T)u*_^(r)«:^(r)+2perms]+c.c., (46) 



— )• — )• 



In all these contributions above, we can substitute the explicit forms of 4^u% into the equation above to get a lot 

ar ^ 

of integrals with r. It's straightforward but the result is rather boring and page-wasting, so we would put them into 
Appendix B. Actually, it has 9 terms differing from each an order of r sequently, plus permutations and complex 
conjugates. At the end we will compile all these terms coming from contributions of all the terms in the Hamiltonian 
according to their indices in a clearer form in order to make our study comfortable. 



D. the shapes of bispectrum 



From Eqs. (|43ll47p and also (jBlllBSp . we can integrate them out to have different shapes. Since the power-law 
indices of t in each integration depend on the value of /i, one may worry that for indices less than —1, the infrared 
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(IR) divergence will occur. However, that is not the case. As has already been shown in Sec. Ill B, the mode solution 
is frozen outside of horizon and there will be no more IR evolution. Actually, the IR divergences will all be canceled 
with each other and we do not see any real singularity, though it is tedious to check analytically and even difficult 
numerically. The same argument can be found in the paper JCAP 1004, 027 (2010) as cited in [29*1. [62*1 Furthermore, 
we can replace the value of u-^ at far future r = to that at horizon-crossing time. Taking these into consideration, 
we first write down all the possible shapes of the bispectrum. Here we define 

< Kir, fct)C(T, fc^)C(T, fc^)| >= (27r)353(^ t)B{h,k2, (48) 

i 

where B{ki, k2, k^) are the shapes of non-Gaussianity. There are totally 10 shapes at the leading order: 



(27r)5Si7^ cos{3fin)T{3^i - 2) 



82/|o(-l)c''' (fclfe^3)'^+^if3M-2|e^ + ^ - l|3r3(-^ + 1) 



(l + ^)(3^2^Al!(l±il)_12(;,_ 1)2) ^^2^2^(1 + ^) , (49) 

^ i>j ^ i 



{2tt) 5 Si72 cos{3fin)r{3fi - 1) 



82/3o(-l)c^^"'(fclfc2fc3)'^+li^3M-l|e, + ^ - l|3r3(-;i + 1) 



l. M'(l + /i) 

2c2 



- Q(P m (n y) E(^'^'I - ^-^^l) + -(6/^' + ^^^^ i2(M - 1)^) y: fc^fc. ]m 



^3.-. = 2(2.)ii/^cos(3M.)r(3M) 4 ( n r 1 (E ^^^1 - E - + 2A)(kl) 

82/3„(-l)cf''-«(fc,fc,fc3).A'3M|,, + /,„l|3r3(-^+ 1) (^2c2 [i-} k, ) ^2^^ ^ ^ 2-^ ^ > ^ ^ 

_ (27r)lM^(M - + M)^g' cos(2M^)r(^ - 2) ^ 

83c^-Vlo(-l)(n. + M - l|^r3(-;, + i) V ^ ' ^ ^ 



^ 2(27r)iM4(^ _ ^)(^ ^ ^)2g4 cos(2A.7r)r(^ - 1) ^ , 

83cf-v^„(-i)(a )i^-^k. + M - +\)^ ^ ^ 



(27r)i(^2 _ i)H^cos{2nn)T{n) 



%^fl,{-l)c''/-^(k,k2k:,Y+^K^\en + ^l - l\^V^{-fi + \) 

,3 



9(M - 1)(2(m - 1)2 - ^2) ^ fc2^2 _ - 1)3 5] fcf + I 5^ hk, I 5^ fcH , (54) 



2c2 

1>J I » \«>J 



2(27r)iij4 cos(2/x7r)r(/x + 1) 



" ^^fl^{-l)cf-\k,k2k^Y+^K^^+^eh + ^l-l\^T^{-^i + \) 



[ (M - 1) f n - i)(2a^ - 1) e ^^^^ + Y E fc') + ^^^v^ E(^'^' - ^»^') 

\ \ i / ^ i>3 i ^ i^J ) 



(55) 



g ^ 18(27r) 2 (^ - l)2g4 cos(2M^)r(A. + 2) 

83/2„(-l)c!^-4(fciA:2A:3)^i^^+2|e,+M-iPr3(-^+i) ' 
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^5.-5 = — (2y^-(l + M)cos(4M^)r(5M-4) ^ ^ _ ^ . 

82/4p(-l)c?^-«(fclfc2fc3)^+^if5^-lQ.+/i-iPr3(-^+i) 1^ ^ ' ' Y ' 

For the next step, we have to discard the terms that wih be divergent as the indices in the original integral becomes 
less than —1, and sum up all the convergent terms to give the total shape. It depends on the value of /i, obviously. 
In the following we give the total shape for different values of /i. The final result is 

Stota/ = ^e(y+l)Sj, , (59) 
y 

where Q is the Heaviside step function defined as: 

{0 , for X < , 
(60) 
1 , for X > 0, 

and y for all the subscripts of B above. From this result, we can also divide all the shapes to four classes: 

i) 'Ssai-I: 'S^-ii '8^1 iB^i+i: Since we assume that the slow-roll parameters eh and ee larger than 0, thus the parameter 
^ > 0, and the indices of these shapes are larger than —1. So they will definitely contribute to Btotai] 

ii) yB3^_3, ;B5p_4, B^^-z'- If /i > 1, i.e., the power spectrum gets a red index, the indices of these shapes are 
larger than —1 and they will contribute to Btotai] 

iii) B^-^: If /i < 1, i.e., the power spectrum gets a blue index, the indices of these shapes are smaller than —1 and 
they will not contribute to Btotai] and 

iv) Bfj_-2- The indices of these shapes are at "divide" values. So when > 1, they will contribute to Btotai while when 
/X < 1, they will not. 

We can also define the estimator through 

B{ki,k2, h) = TFNL{^^Vi^ri^ + 2perms.} , (61) 

so, 

B{ki,k2,k3) 



Fm l ~ — 



Bik^MM (62) 



for each shape listed above. The result is rather obvious by just substituting each shape into Eq. ([5]), so we will 
not list them here in order to save the page. We will only show the equilateral limits of these i^jvL's, of which 
ki — k2 = k^ = k, m. the next paragraph. 



E. The Equilateral Limit (fci = k2 — ki — k) 

In this section, we can take the equilateral limit, namely, ki = k2 — k^ — fc, which gives simpler form of F^l- 
Following the last paragraph, we have: 

{F.l)Z% = fc2-2. 5(2.r)ls]cos(3M.)r(3;.- 2)r(-/.+ l) ^^^(1+^ _ _ 

(F..)-- = 5(2.)iEcos(3,.)r(3, l)r(-,+ l) _ _ 

\ NL)3^.2 24 x 33M-i7r2/™(-l)A^3cs^"'i?2|,,^ + ^_l|^ d W y ) 
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(p -.egml _ ,2-2u 5(27r)2 COs(3/i^)r(3/i)r(-/i + |) 1 

^^^^^^-^ 192x3M-2^2|,^ + ^_l| ' (66) 

,p ,e,u,i 5(27r)i/z^(A. - 1)(1 + fif cos(2M7r)r(M - m-f^ + j) 

(^^^^'^-^ ^ 96x3M-%2|e. + M-l| ' ^^^^ 

^e,u»^ _ 5(27r)i(/x^ - l)cg cos(2M^)r(A.)r(-/i + 1) n„2 , M% 

^^^^^/^-i 8x3^^Vk/. + M-l| ^ ~ ' ' ' ' ^ ' 

(F ^e...^ _ 5(2^)^cg cos(2/i7r)r(/i + l)r(-M + ^) 3 

(^^^l), 24x3M+i^V|e. + /^-l| -18(M-l)(2/i-l)), (69) 

,p ,e,™; 5(2^)1 (/i - l)^c^ cos(2A.7r)r(M + 2)r(~;. + 1) 

(^^^^''+1 32x3A'+Wk.+A^-l| ' ^ ^ 



^ 5(27r)is2(l + fi) cos(4/x7r)r(5Ai - 4)r(-/Lt + i) 
24 X 35A'-4,r2/Io(-l)A''^'c^"6|e;, + ^ - 1| ' 



,e,mi _ ,A-4, 5(2^)^S^ cosiA^^7^)Ti5^, - 3)r(-A. + 1) 



(-f'ArL)5M-5 = « ... .9 . -.^ 4»-6i ^ — ' i'l) 

^ """^^^-^^ 12x35A<-V/lo(-lV^^^c!'^-'k/« + M-l| ' 

and same as the shape, we have for total non-linear parameter: 

(^a^l)XI = E ©(J^ + l)(^^i^)r^' ■ (73) 

V 

From the result above we can see that (FnlT''"'' has a slight running behavior due to the deviation of fj, from 1, 
that is, due to the non-minimal coupling behavior. This is different from the usual minimal coupling case which was 
studied in 27] where the equilateral limit of Fnl is independent of k. This is another result in this paper and will 
later be confirmed with numerical calculations. One can also obtain F^l in the local limit (fci « fc2 ^ ^3) and folded 
limit {ki = 2/c2 — "^k^). Whichever limits they are in, all the values of F^l with different indices can also be divided 
into four classes by the same criteria used for B{ki, ^2, fca)- 

F. An explicit example: non-minimal coupled chaotic inflation 

In order to support our long analytical derivation, we in this section focus on an explicit model of non-minimal 
coupling inflation. For simplicity but without losing generality, we consider chaotic inflation, of which the potential 
has a quadratic form as V{(j)) = A(/)^/4 where A is the coupling coefficients. Furthermore, we set the non-minimal 
coupling term f{R, 4>) = + S,R4>'^ where G is the Newtonian gravitational constant and ^ is the non-minimal 
coupling coefficient. This model is indeed very interesting since with the presence of non-minimal coupling term, the 
coefficient A doesn't need to go to incredibly small value (Fakir and Unruh, 0) to meet the observational constraint, 
and this property has been used to construct Higgs inflation models, which connects inflation theory to particle 
physics in the Standard Model [52]. In Higgs inflation, the Higgs potential asymptotically coincides with the chaotic 
potential in inflation period where the scalar field is in a high energy region with some large value, and a too small 
value of A will make the models inconsistent with the constraints from Standard Model [Hj . Putting aside its physical 
motivation, in this paper we study whether it can give rise to large non-Gaussianities. 
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From the original action ([Ij, one can obtain the equation of motion for the fields: 

4> + 3H^-6aH + 2H^)cj,+ ^^^ = , (74) 

and Friedmann equation; 

+ icj,^) - + V{<t>) ~ e^Hcjy^ , (75) 
and the above two equations can be combined to give another equation: 

H = + e0') = -l'}>^ + iH^<p ~ ~ m ■ (76) 

Firstly, we draw the background evolution of the system in Figs. [1] and [21 From the plots we can see that with 
natural choice of initial conditions and parameters, a period of inflation can be easily constructed with enough amount 
of number of e-folds. Setting ^ = 1000, the parameter A could be raised up to 0(10"^), compared to the unnatural 
choice of A ~ 10~^^ in the case with ^ = (Fakir and Unruh, Q)- Figs. [3] and |4] shows the amplitude and the 
fc— dependence of its quadratic perturbation spectrum, whose analytical form has already been given in Eq. p7p . For 
the given initial conditions and parameters, we can see that the spectrum behaves nearly fc- independent, with a slight 
tilt caused by the deviation of \x from 1, which is mildly favored by the WMAP-7 data [53|. The amplitude of the 
spectrum is also consistent with the observations. 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 
-0.2 




FIG. 1: The evolution of slow roll parameter e w.r.t. cosmic time t. The arrival of e at 1 stops the inflation. Parameters and 
initial values: ^ = 1000, A = 10~^, 0i = 4.9, 0, = 0.063. The normalization is SttG = 1. 

Next let's move on to the non-Gaussianities that this model can give rise to. Since the power spectrum of this model 
has a blue tilt, the first two classes of the total four in the shape B{ki, k2, k^) as well as the estimator Fjvl which 
was shown in the last paragraphs will be applied. From our numerical calculations, we can obtain the values of every 
parameter that appears in Eqs. P^ -([55 | as well as (l551) - ([7^ . With this in hand, we can easily obtain the numerical 
results of the non-Gaussianities generated in this model. The total shape (of leading order) of the non-Gaussianities 
and the estimator in the equilateral limit {FNLy™^ are shown in Figs. [S] and [SI From the plots we can see that 
the shape of the non-Gaussianities are well within the constraints of the observational data. Note that the estimator 
(-Fjvl)^*"*' shows a running behavior with respect to fc(= fci = fc2 = fcs), with a positive sign. This is because of the 
effect of non-minimal coupling which makes the parameter ^ deviate from 1 and thus {FmlY''^^^ will be dependent 
on k. In the GR limit /i ^ 1, {FnlY'^^^'' will have a constant value, as shown in [27| . 



IV. CONCLUSION AND DISCUSSIONS 



In this paper, we performed the non-Gaussianities of a general single scalar field which linearly couples to gravity. 
Our result shows that due to the non-minimal coupling, the power spectrum will deviate from scale-invariance, which 
in order lead to the complicated non-Gaussianities in the 3rd order. We obtained all the possible shapes of the 3-point 
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FIG. 2: The e-folding number A'^ w.r.t. cosmic time t. Parameters and initial values are the same as in Fig. [T] 



|p| 1E-11 




FIG. 3: The amplitude of power spectrum w.r.t. cosmic time t. 




FIG. 4: fi w.r.t. cosmic time t. Slightly deviation from 1 is obtained due to the non-minimal coupling. 
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FIG. 5: The shape of non-Gaussianities B{k\, k2, fea). Here we renormahze x = ki/ks, y = fc2/fc3, and set fca — 1. 



F_NL 



37 - 




FIG. 6: The estimator of non-Gaussianities in the equilateral limit; (Fjvl)'^''"''- Running behavior is obtained due to the 
non-minimal coupling effect. 

correlation functions and for different tilt of power spectrum, we showed that different shapes will be involved in to 
give rise to non-Gaussianities. Our calculation presents the description in general non-minimal coupling inflation and 
this result, if verified to all the orders, can provide relation between 2- and 3-point correlation functions and can be 
used to constrain non-minimal coupling models. 

Another result that was presented in this paper is that there is some running behavior of the estimator Fn l in the 
equilateral limit with respect to k which is different from the normal minimal coupling case. This behavior is due to 
the non-minimal coupling, and are expected to have signature on observations in order to distinguish minimal and 
non-minimal cases. 

Besides the analytical calculations, we also performed numerical computations on a specific example of non-minimal 
coupling chaotic inflation model. This model is extendedly studied with the application of Higgs inflation. We obtained 
the behavior of background, 2-point power spectrum as well as the shape and estimator of its non-Gaussianity. We 
showed that the non-Gaussianities are well within the observational constraint, with the running behavior of F^^^'' 
w.r.t. k. 

Other than inflation, such non-minimal coupling system can also be applied to other aspect in cosmology. For 
example, non-minimal coupling theory can act as dark energy [55| or give rise to a bouncing/cyclic universe (56|. 
Moreover, non-minimal coupling can be used to make up an open/closed universe [57| . while dualities of Einstein's 
gravity in the presence of a non-minimal coupling was taken on in [58] . The stabilities and singularities in superac- 
celeration phases was discussed in [s^ , and the removal of singularities in Loop Quantum Gravity with non-minimal 
coupling was studied in j6Q||. Our calculation of non-Gaussianities are also expected to be applied to these interesting 
fields. 
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Appendix A: extend to non-linear coupling 

Here, we only list the constraint equations for the general form of non-minimal coupling single scalar field, including 
the non-linear coupling case. The constraints from N and Nt are: 

4aE + l2HfRoC - UH^fBoa - Aa-^fnoHd^^ 

~Aa-^fBod\ - [UaH - 6C + 2a-^d^i;)dtfm + 72{H - H^){H + 2H^)fRBoa - 36(ij - H^)fRmQ 
-36(ij - H^)HfRm{K + I2a-^{H - H^)fRRo{d'^^ + 2Hd^i^ + d^a + 2d\) 
+6Hdt[-l2{H + 2H^)fRRoa + 6HfRRo{4C - a) + GfRRoC - 2a-^fRRoid^ + 2Hd''^lj + d'^a + 25^)] 
+24a-\H + 2H^)fRRod^a - Ua'^ H fRRod^AC -a)- I2a~^ /rroO'^C 

+^a-'^fRRod^{d^ip + 2Hd^i} + d^a + 23^ Q = (Al) 



-^a-^fRod'^N' - \2a-^H{H + 2H^)fRRod'a + Qa'^ f RRad\K - a) + Qa'^ H f RR^d' '( 

~2a-^HfRRod\d^ + 2i/aV + d^a + 2d^0 + a-^dt[12{H + 2H^)fRRod'a ~ 6HfRR,od\4C - a) 
-QfRBodX + 2a-^fRRod\d^ + 2Hd^^j + d^a + 2d\)] 

+2a-^fRoHd'a + a-^dtfRod'a - 2a~^fRad\ ^ (A2) 

respectively, where Jrq and fRRa denotes the background value of the first and second derivatives of /(i?, (j)) with 
respect to R. 



Appendix B: Contributions from terms in //f„j w.r.t. r 

Here, we list the contributions from terms in Hf^^, i.e. Eq. (|43ll47p . in terms of r by substituting jifU^ in. One can 
see that there contains integrals of different power-laws of r, each differing one order from the other in every term. 
One can combine the integrals of the same order power-law to have more neat forms, as in Sec. Ill D. 

The contribution from C,^: 

< C3 > ^ g (2.)3^3(^^fc>)(s + 2A).lg^e-3-^ 1^ d.e--nd3M- 

^"^pR^{~l)cs''~\k^k2k:iY\eh + ^l-\?T^{^^l + \)J-oo 

8i/2„(-l)c^-6(fcifc2fc3)^|eh + A* - iPr3(-^ -t- 1) 
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The contribution from C^'^: 



8i /3o(-l)c?''(/eiA;2A;3)''+2|e^ + /X - l|3r3(-/x + i) 
X / rfTe''=»'^^[{+4c>(l + /i)fc|fc^|T|3''-3+8icffcifcifc^|T|3''-2} + 2perms.] 

J —OO 
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8ip^,{-i)cl''{hk2hY+^\eh + i|3r3(-^ + 1) 

X / (iTe''=-^^[{+4cXl+At)A:|A;||T|^-i+8ic^fciA;|fc||T|''} + 2j3erms.]+c.c. (B2) 

J —OO 

The contribution from C(^C)^- 

+2ic,M^(l + /x)2K|r|'^-2 - 4clfx{l +IJ')\Y, ^i^i ) 1^1""' " 8ic3(fcifc2A:3)|r|'^} + 2perms.] + c(<fi3) 



The contribution from 

I ■>i'>'rrV>A"l N 



2^ 12(2 ^)3^3(^^ fc07ri(A^- l)2g4e-3.iV 

3^-2 



J —00 



-00 

4(27r)353(^. fct)7ril]if2g-3i^ 



8^/io(-l)c«''"'(fcifefc3)''+2|e/. + M - l|3r3(-/, + i) 

/ rfrp2 • A;3)fc?e*'=^-^^[M2(i ^ /x)2|r|3/'-3 + 2ic,M(l + m)(A;2 + fc3)|Tp''-2 - 402^2^3 Irp""!] + 2perms\ 

J — 00 



X 

-00 

+C.C. (B4) 



The contribution from C,dC,dx'- 



<CdCdx> 3 . ..^3, 



18(27r)353(^^ ki)nHfi-iyH^e- 



3i 



^/™(-i)cr(fcifc2fc3)^+V|e/. + M-ipr3(-/i + i) 

/ dTe^'=^^^[^^2_^i{+4cXl + ix)kfki\T\^-^ + 8iclk2k''^klm + 5perms] 
J-00 '^3 



IRO 
rO 

X 

-co '''3 
^ 2(27r)3^3(^,fc>)^|S2g-3»iy 

85/io(-l)c«''(fcifefc3)''+V|eh + M - l|3r3(-/, + i) 

X /° dre*'^-^-[.^H_^{+4e4^,(l + ^)fc2fc2|^|5p-5 ^ 8ic^A:2fc?A:||T|5'^-4} + 5perms] 
J-00 ^3 

12(2^)353(^^ fct)^i(M - l)2Si/2e-3«^ 



8f/3o(-l)cf^(fcifc2fc3)''+V|e/. - l|3r3(-/x+ i) 



/° dTe''''^^['^^-p-{+Acin{l + fj.)klkl\Tf^-^ + 8ic%klkl\Tf^'-^} + 5perms] + c.c. (B5) 

J-oc 1^3 
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